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Abstract 

The resonant substructure in D° —> Kgir + ir~ decays is described by a combination of ten quasi 
two-body intermediate states which include both CP-even and CP-odd eigenstates and one doubly- 
Cabibbo suppressed channel. We present a formalism that connects the variation in D° decay time 
over the Dalitz plot with the mixing parameters, x and y, that describe off-shell and on-shell D°—D° 
mixing. We analyze the CLEO II.V data sample and find the parameters x and y are consistent 
with zero. We limit (—4.7<x<8.6)% and (—6.1 <y<3.5)% at the 95% confidence level. 

PACS numbers: 13.25.Ft, 12.15.Mm, 11.30.Er, 14.40.Lb 
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Studies of the evolution of a K° or B° into the respective anti-particle, a K° or B °, 
have guided the form and content of the Standard Model and permitted useful estimates of 
the masses of the charm [1] and top quark [2| prior to their direct observation. A D° can 
evolve into a D° through on-shell intermediate states, such as K + K~ with mass, uik+k- — 
moo , or through off-shell intermediate states, such as those that might be present due 
to new physics. This evolution through the former (latter) states is parametrized by the 
dimensionless variables —iy (x) defined in Eq. [221 

Many predictions for x in the D° —> D° amplitude have been made Q. Several non¬ 
standard Models predict |x| > 0.01. Contributions to x at this level could result from the 
presence of new particles with masses as high as 100-1000 TeV [ 5 J. The Standard Model 
short-distance contribution to x is determined by the box diagram in which two virtual 
quarks and two virtual W bosons are exchanged. The magnitude of x is determined by the 
mass and Cabibbo-Kobayashi-Maskawa (CKM) [§] couplings of the virtual quarks. From the 
Wolfenstein parameterization |7] where A = sin 2 9q ~ 0.05, contributions involving b quarks 
(~A 6 ) can be neglected relative to those with d and s quarks (~A 2 ). The most prominent 
remaining amplitude is proportional to (m 2 — The near degeneracy on the W 

mass scale of the d and s quarks results in a particularly effective suppression by the GIM@ 
mechanism. A simple estimate of x is obtained by comparing with the Kaon sector; 

AM D o _ f D o (m 2 - m 2 d )m D 0 
AMft'O f K o{rn 2 c - m 2 u )m K 0 ' 


Assuming f D 0 ~ f K 0 and taking m u = 5 MeV, rrid = 9 MeV, m s = 60 — 170 MeV, 
m c = 1.2 GeV and AMj-o = (3.48±0.01) x 1CT 15 GeV, and x = Afi ^ D ° = 6.31 x 10 11 x AMpo 
yields, x = 2 x 10 _5 -2 x 10~ 4 . Short distance contributions to y are expected to be less than 
x. Both are beyond current experimental sensitivity. Long distance effects are expected to 
be larger but are difficult to estimate due to the large number of resonances near the D° 
pole. It is likely that x and y contribute similarly to mixing in the Standard Model. Decisive 
signatures of new physics include \y\ -C |x| or Type II or Type III CP violation [9]. I 11 order 
to assess the origin of a D°—D° mixing signal, the values of both x and y must be measured. 
Previous attempts to measure x and y include: the measurement of the wrong sign 

which is sensitive to the mixing rate Rm = 
K + K~, 7r + 7r _ [l2| which are sensitive to y\ 
hadronic branching ratio which measures 


semileptonic branching ratio D° —> K£u 
decay rates to CP eigenstates D° 

K+ 7T 


P+y 2 


and the wrong sign D° 


la, M, 115 


„/2 


x'“ = (r/sinh^Tr + ^coshxTr) 2 and y' = y cos 5 k-k — % sin Sk-k- Here, dxn, which has yet to 
be measured experimentally, is the relative strong phase between D° and D° to K + tt~. 
In this study we utilize the fact that the values of x and y can also be determined from 
the distribution of the D° —> Kg7T + n~ Dalitz plot if one measures that distribution as a 
function of the D° decay time. We show that x and y can be separately detemined. This 
is the first demonstration of possible sensitivity to the sign of x. Predictions of the sign of 
x are sensitive to the details of the treatment of long distance effects within the Standard 
Model as well as the nature of potential new physics contributions. 

The time evolution of the D°-D° system is described by the Schrodinger equation 


dt\D°{t)J 



( 2 ) 


where the M and T matrices are Hermitian, and CPT invariance requires M n = M 2 2 = M 
and Tn = T 22 = r. The off-diagonal elements of these matrices describe the dispersive or 
long-distance and absorptive or short-distance contributions to D°-D° mixing. 
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The two eigenstates D\ and D 2 of the effective Hamiltonian matrix (M — |T) are given 


by 

\D 1 , 2 )=p\D 0 )±q\D 0 ) , p 2 + q 2 = 1. 
The corresponding eigenvalues are 

Ai,2 = ^i,2 — 2 44,2 = -T^j ^M 12 — -T 12 

where mi j2 , T 12 are the masses and decay widths and 


Q 

P 


M { 2 - A 


2 1 12 


\ Mn - |T 12 ' 

The proper time evolution of the eigenstates of Eq. [2] is 

\D\ )2 (t)) = eip(t)\Di t2 )i e\ j2 (f) = e 
A state that is prepared as a flavor eigenstate |D°) or \D°) at t — 0 will evolve according to 




(3) 

(4) 

(5) 

( 6 ) 


D °(t)) = 7 ^ [p( e i(t) +e 2 (t )) | D°) +q(e{t)-e{t))\D°) 

I D°(t))= ^[p(e{t)-e{t))\D°)+q(e{t)Ait))\D°) 


(7) 

( 8 ) 


We parameterize the Dalitz plot following the methodology described in Ref. E 


17] using the same sign convention as Ref. [181. Il9l. 20]. Now, however, we generalize to the 


case where the time-dependent state is a mixture of D° and D° so the Dalitz Plot distribution 
depends also on x and y. We express the amplitude for D° to decay via the j'-th quasi-two- 
body state as aje lSj A 3 k where A], = A J k (m 2 K07r -,rri 2 n ) is the Breit-Wigner amplitude for 
resonance j with spin k described in Ref. [T3] - We denote the CP conjugate amplitudes for 
D° as A{ = Al(m 2 K o n+: m^). 

We begin our search for D° — D° mixing in D° —> Kg7r + 7r~ from the results of our 
standard fit in Ref. E] which clearly observed the ten modes, ( K*~n + , Jlg(1430) _ 7r + , 
Jl*(1430)-7T + , /W(1680)-7 t+, K° s p, K° s cu , K° s f 0 ( 980), K° s f 2 ( 1270), K° s f 0 ( 1370), and the 
“wrong sign” K* + tt~) plus a small non-resonant component. 

The decay rate to KgTT + n~ with [yn 2 K o^_ 1 m 2 +7r _) at time t of a particle tagged as |D°) 
at t — 0 is 


dr(m 2 K 7 r ,ml n ,t) 


■\M\ 2 dm 2 Kn dm: 


2 

7T7T "> 


256vr 3 Af 3l " N ^ 

where the matrix element is defined as A4 = (f\7i\i). We evaluate |A4| 2 where | i) is given 
by Eq.H and (f\ = {K° s ir + Tr-(m 2 K o^, m 2 +7 r _)l- 

The decay channels can be collected into those which are CP-even or CP-odd (with 
amplitudes A + or A_) and to those which are D° or D° flavor eigenstates (with amplitudes 
A f or A f y, 


</|W|£> + ,_) 

— ^2 a j e J ‘^fc=even,odd ~ ^+,- 

(10) 


= a t e ^' 7 *7li=eveii,odd = A+,- 

(ii) 

m\D° f ) 

= -4/ 

(12) 

mm) 

= =a t 

(13) 
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Dalitz plot analyses are sensitive only to relative phases and amplitudes. As in Ref. 
we fix a p — 1, S p — 0 and assume a 3 = dj, Sj = Sj. In Ref. [20|, we considered CP violation 
more generally and allowed a 3 ^ dj, Sj ^ Sj. 

Collecting terms with similar time dependence we find 


(f\H\D°(t)) = i((/|W|0,M) + {/|W|£> 2 («))) (14) 

= ^((/|W|(pC° + ?I5 0 )) e iW 

+ (/|W|(pD‘>- g D»)»e 2 (() 

= ^([p( A f+ A + + A_)+q(Af+A + +A_)]ei(t) 

+ \p( A f + A + + A -)~q( A f + A + + A_)]e 2 (t)) 

= 2 ^ + Xf)Af+(l + X+)A++(1 + X-) A -] e i(t) 

+ 2 ^ _ Xf)Af+(l ~ X+)^++(l - X-) A -] e 2(t) 

= ei(t)Ai + e 2 (t)A2 

(j\nm*))= A(mmt)) - (mmtm (is> 

2 q 

= 2[(^+Xj 1 )Aj+(l+x+)A + -{l+X-)A-\ e 1 (t) 

+ l[(i-X7 1 %+(i-x: 1 )3--(i-x: 1 M-]e 2 (i) 

= £i(t)A\ + e 2 (t)A 2 , 

for D° and D°, respectively. Similar to Ref. [sj, 


Xf 


*7 


x± 


q A f — 

A f 

1 — e 

p A f 

A f 

l + e 

q a 7 

A 1 

1 — e 

pAj 

A f 

l + e 


= ±I— i e W 

p a ± l + e 


e i{5+<t>) 


e -i(6-<t>) 


(16) 

(17) 

(18) 


where S is the relative strong phase between D° and D° to Kg7i + n~, and in the limit of CP 
conservation, the real (7P-violating parameters, e and 0, are zero. Squaring the amplitude 
and factoring out the time dependence yields 


\M\ 2 = \ei(t)\ 2 \Ai\ 2 + \e 2 (t)\ 2 \A 2 \ 2 


M 


+2^[e 1 (t)e*(t)A 1 A* 2 \ 
— |ei(f)| 2 A\ " + |e 2 (7)| 
+23? [ei (t) (t) A\ A 2 \. 


Ar 


(19) 

( 20 ) 


The time-dependent terms are given explicitly by 

l e i ,2 (7) | J = exp(2A(Ai >2 )t)=exp(-r 1>2 f) (21) 
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ei(t)e 2 (t) 


( 22 ) 


= exp (—T(1 ± y)t) 

= exp (— iXit) exp ( +i\ 2 t ) 
= exp (— T(l+ix)t) 


where 


r = 


ri+r 2 


X = 


771 1 — m2 


y = 


ri-r 2 


2 r ’ " 2r ' ^ 

Experimentally, y modifies the lifetime of certain contributions to the Dalitz plot while x 
introduces a sinusoidal rate variation 

This analysis uses an integrated luminosity of 9.0 fb _1 of e + e~ collisions at A /s ~ 10 GeV 
provided by the Cornell Electron Storage Ring (CESR). The data were taken with the 
CLEO II.V detector 
of D° 


K° s tt h 


7T 
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The event selection is identical to that used in our previous study 
IS, 20] which did not consider D°—D° mixing. We reconstruct candidates 


TTcj or 


for the decay sequence D* + —> D u —> Kg7i + 7i~. The charge of the slow pion 

7 Tg) identifies the initial charm state as either D° or D°. The detector resolution in the 
Dalitz plot parameters m 2 Kn and m 2 ^ is small relative to the intrinsic widths of intermediate 
resonances; the exception is the decay channel D° —> KgU), oj —> 7r + 7r _ . We reconstruct the 
D° decay time t as described in Ref. |l3| ] . 

The uncertainty in t, a t , is typically 200 fs or 0.5/T and cannot be neglected. We fit 
the unbinned decay time distribution by analytically convolving the exponentials in each 
term in Eqs. [201 and ED by a resolution function similar to, but slightly modified from, 
that used in Ref. 0 and Ref. [2lJ. The signal likelihood is represented as the sum of an 
exponential convolved with two Gaussians. The width of the first Gaussian is the event-by- 
event measured proper time error, cr t , times a scale factor, S$ i g , which allows for a uniform 
mistake in the covariance matrix elements of the D meson and its daughters, perhaps due to 
an imperfect material description of the detector during track fitting. For the other Gaussian, 
the measured proper time errors are ignored and the width OMisSig and the normalization 
/MisSig are fit for directly. This Gaussian models the ‘MIS-measured SIGnal” proper time 
resolution when the measured a t is not correct, as would be the case for hard multiple 
scattering of one or more of the D meson daughters. The sum of these two components to 
the likelihood is normalized by the total signal fraction / sig . Note that if we understand our 
detector well, we will find that the scale factor used in the first Gaussian is close to unity 
and the fraction of the signal in the second Gaussian is near zero. 

The treatment of the background is similar to that of the signal. The total background 
likelihood is normalized by the background fraction, which is (1 — / S i g ). We consider two 
types of background: background with zero lifetime and background with non-zero life¬ 
time tbg normalized by / TBG . We constrain both backgrounds to have the same resolution 
function. The model for the resolution function is two Gaussians, with core width Obg, 
misreconstructed width (TMisBG and the background fraction /misBG bi the wider Gaussian. 

We perform an unbinned maximum likelihood fit to the Dalitz plot which minimizes the 
function T given below 

T = Y,~ 21n£ + ^-21nZ, (24) 

D° TJo 


where C and C are defined as in Ref. 20] using M. and M. as defined in Eqs. [20] and 
convolved with the resolution function described above. Simplified Monte Carlo studies 
indicate that our fit procedure is unbiased and the statistical errors as determined by the fit 
are accurate. 
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FIG. 1: Projection of the results of Fit A onto the D° decay time with a) linear and b) logarithmic 
vertical scale. 


Our standard fit to the data, described above, is referred to as Fit A. Fit B is identical 
to Fit A except CP conservation (e = 0, (j) = 0) is assumed. The D° and D° sub-samples 
are fit independently in Fit Cl and Fit C2, respectively. Fit Cl and Fit C2 are identical to 
Fit B. 

Fit A has 35 free parameters; ten resonances and the non-resonant contribution cor¬ 
respond to ten relative amplitudes and ten relative phases, signal fraction and mis-tag 
fraction, four signal decay time parameters, five background decay time parameters, two 
mixing parameters and two CP-violating parameters. The results for x, y, e and </> are 
in Table U and are consistent with the absence of both D° — D° mixing and CP violation. 
The one-dimensional, 95% confidence intervals are determined by an increase in negative 
log likelihood (—21n£) of 3.84 units. All other fit variables are allowed to vary to distinct, 
best-fit values. The amplitude and phase, a,j and Sj, for all fits in Table Q] are consistent 
with our “no mixing” result [3]. The projection of the results of Fit A onto the D° decay 
time is shown in Fig[0 

We find the parameters describing the signal decay time, / sig = (97.1 ± 0.8)%, r sig = 
402 ± 6 fs, S sig = 1.13 ± 0.02, (x M isSi g = 735 ± 155 fs, (1 - f MisSig ) = (96.9 ± 1.5)% and 
the parameters describing the background time, / Tbg = (100 ± 1)%, Tb g = 94 ± 59 fs, 
(1 — /misBg) = (87 ± 11)%, o ’bg = 197 ± 39 fs, o’MisBG = 1116 ± 321 fs. The scale factor 
S sig, although not consistent with unity, is comparable to results from other CLEO lifetime 


analyses which include Ref. [Ill, [13|, 21 


We evaluate a contour in the two-dimensional plane of y versus x that contains the true 
value of x and y at 95% confidence level (C.L.) without assumption regarding the relative 
strong phase between D° and D° Kg7r + 7i~. We determine the contour around our best-fit 
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TABLE I: Results of the Dalitz-plot vs decay time fit of the D° —> KgTT + TT~. Fit A allows both 
D° — D° mixing and CP violation. Fit B is the CP-conserving fit, e = 0 and 0 = 0. Fit Cl 
(C2) is the fit to the D° ( D °) sub-sample. The errors shown for Fit A and Fit B are statistical, 
experimental systematic and modeling systematic respectively and the 95% confidence intervals 
include systematic uncertainty. The errors for Fit Cl and Fit C2 are statistical only. 


Parameter Best Fit 1-Dimensional 95% C.L. 


Fit A 

Most General Fit 

x (%) 

2.6ljj;g±0.4±0.4 

|x| < 9.8% 

y(%) 

—O.attie ± 0-8 ± 0.4 

| y\ < 9.5% 

e 

-0.3 ±0.5 


0 (°) 

42 ±78 


Fit B 

CP-conserving fit 

x (%) 

1.9j;jj;§ ±0.4 ±0.4 

(-4.7:8.6) 

y(%) 

-1.4 ±2.4 ±0.8 ±0.4 

(-6.1:3.5) 

Fit Cl 

D° sub-sample 

x (%) 

o q+5.0 

O.O— 4 8 

(-6.1:13.5) 

y(%) 

q o-|-3.6 

-z.z S7 

(-10.2:4.2) 

Fit C2 

D° sub-sample 

x (%) 

0-6±|;I 

(-16.0:11.5) 

y(%) 

-o.tti 

(-6.6:13.0) 


values where the —21n£ has increased by 5.99 units. All fit variables other than x and y 
are allowed to vary to distinct, best-fit values at each point on the contour. The contour 
for Fit A is shown in Fig. [21 On the axes of x and y, these contours fall slightly outside 
the one-dimensional intervals listed in Table Q] as expected. The maximum excursion of the 
contour of Fit A (Fit B) from the origin corresponds to a 95% C.L. limit on the mixing rate 
of R m < 0.84% (R m < 0.55%). 

We take the sample variance of x, y, e and 0 from the nominal result compared to the 
results in the series of fits described below as a measure of the experimental systematic and 
modeling systematic uncertainty. 

We consider systematic uncertainties from experimental sources and from the decay model 
separately. Our general procedure is to change some aspect of our fit and interpret the 
change in the values of the mixing and CP-violating parameters in the non-standard fit 
relative to our nominal fit as an estimate of the systematic uncertainty. Contributions 
to the experimental systematic uncertainties arise from our model of the background, the 
efficiency, the event selection criteria, and biases due to experimental resolution as described 
in Ref. [19j. Additionally, we vary aspects of the decay time parametrization. To estimate the 
systematic uncertainty regarding the uu, dd, ss content of the background, we perform fits 
where the background is forced to be all zero lifetime and all non-zero lifetime. We consider 
a single or a triple rather than a double Gaussian to model the decay time resolution of the 
signal and background. We also vary by ±lcr the fraction of misreconstructed signal /MisSig- 
Finally, we set the scale factor for the measured proper time errors S sig to unity. Variation in 
the event selection criteria are the largest contribution to the experiment systematic error. 

Contributions to the theoretical systematic uncertainties arise from our choices for the 
decay model for D° —> Kg7r + 7T~ 


as described in Ref. lfij]. We also consider the uncertainty 












arising from our choice of resonances included in the fit. To study the stability of our results 
with other choices of resonances, we performed fits which included additional resonances to 
the ones in our standard fit. We compared the result of our nominal fit to a series of fits 
where each of the resonances, a or / 0 (600) and / 0 (1500) which are (TP-even, and p(1450) and 
p(1700) which are (TP-odd were included one at a time. In the standard fit we enumerate 
the non-resonant component with the AT* resonanaces. We also considered fits where the 
non-resonant component was considered to be (TP-even or (TP-odd. Finally, we consider 
a fit that includes doubly-Cabibbo suppressed contributions from Ao(1430), AT(1430) and 
A'* (1680) constrained to have the same amplitude and phase relative to the corresponding 
Cabibbo favored amplitude as the A'* (892). There is no single dominant contribution to the 
modeling systematic error. 

In conclusion, we have analyzed the time dependence of the three-body decay D° —> 
Kg7i + n~ and exploited the interference between intermediate states to limit the mixing 
parameters x and y without sign or phase ambiguity. Our data are consistent with an 
absence of both D°—D° mixing and (TP violation. The two-dimensional limit in the mixing 
parameters, x versus y, is similar to previous results obtained from the same data sample (l3| . 
when assumptions regarding Sk-k are removed. We limit (—4.7 < x < 8.6)% and (—6.1 <y< 
3.5)%, at the 95% C.L., with the assumption of (TP-conservation. We measure the CP- 
violating parameters e = —0.3 ± 0.5 and (ft = 42° ± 78°. 

We thank Alex Kagan, Yuval Grossman, and Yossi Nir for valuable discussions. We 
gratefully acknowledge the effort of the CESR staff in providing us with excellent luminosity 
and running conditions. This work was supported by the National Science Foundation, the 
U.S. Department of Energy, and the Natural Sciences and Engineering Council of Canada. 
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FIG. 2: Allowed regions in the plane of y versus x. No assumption is made regarding S K o 7T + 7T ~. 
The two-dimensional 95% allowed regions from our Fit B (light shaded region) is shown. The 
allowed region for Ar is the average of the yep 0,0 results. Also shown is the limit from 
D°—*K(*)Iv from BABAR 10]. All results are consistent with the absence of mixing. The limits 
from CLEO 0 ] and BABAR 14] from D —► Kn have similar sensitivity to Fit B. The 95% allowed 
regions (not shown) are circles of radius 5.8% and 5.7%; respectively, when assumptions regarding 
5Kn are removed. The 95% allowed region from Belle [151] also from D —► A"7r is more restrictive - 
a circle of radius 3.0%. 
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